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The field of asteroid thermophysical modeling has experienced an extraordinary growth 
in the last ten years, as new thermal infrared data became available for hundreds of thousands 
of asteroids. The infrared emission of asteroids depends on the body’s size, shape, albedo, 
thermal inertia, roughness and rotational properties. These parameters can therefore be derived 
by thermophysical modeling of infrared data. Thermophysical modeling led to asteroid size 
estimates that were confirmed at the few-percent level by later spacecraft visits. We discuss 
how instrumentation advances now allow mid-infrared interferometric observations as well as 
high-accuracy spectro-photometry, posing their own set of thermal-modeling challenges. We 
present major breakthroughs achieved in studies of the thermal inertia, a sensitive indicator for 
the nature of asteroids soils, allowing us, for instance, to determine the grain size of asteroidal 
regoliths. Thermal inertia also governs non-gravitational effects on asteroid orbits, requiring 
thermophysical modeling for precise asteroid dynamical studies. The radiative heating of 
asteroids, meteoroids, and comets from the Sun also governs the thermal stress in surface 
material; only recently has it been recognized as a significant weathering process. Asteroid 
space missions with thermal infrared instruments are currently undergoing study at all major 
space agencies. This will require a high level of sophistication of thermophysical models in 
order to analyze high-quality spacecraft data. 


1. INTRODUCTION 

Asteroid thermophysical modeling is about calculating 
the temperature of asteroids’ surface and immediate sub¬ 
surface, which depend on absorption of sunlight, multiple 
scattering of reflected and thermally emitted photons, and 
heat conduction. Physical parameters such as albedo (or re¬ 
flectivity), thermal conductivity, heat capacity, emissivity, 
density and roughness, along with the shape (e.g., elevation 
model) of the body, its orientation in space, and its previous 
thermal history are taken into account. From the synthetic 
surface temperatures, thermally emitted fluxes (typically in 
the infrared) can be calculated. Physical properties are con¬ 
strained by fitting model fluxes to observational data. 

One differentiates between sophisticated thermophysi¬ 
cal models {TPMs; Lebofsky and Spencer, 1989; Spencer, 
1990; Spencer et al., 1989; Lagerros, 1997, 1996a, 1998; 
Delbo, 2004; Mueller, 2007; Rozitis and Green, 2011) and 
simple thermal models, which typically assume spherical 


shape, neglect heat conduction (or simplify its treatment), 
and do not treat surface roughness (see Harris and Lagerros, 
2002; Delbo and Harris, 2002, for reviews). In the past, us¬ 
age of TPMs was reserved to the few exceptional asteroids 
for which detailed shape models and high quality thermal 
infrared data existed (Harris and Lagerros, 2002). In the 
last ten years, however, TPMs became significantly more 
applicable (see § 6), thanks both to new spaceborn infrared 
telescopes {Spitzer, WISE and AKARI; see Mainzer et al., 
2015) and to the availability of an ever-growing number of 
asteroid shape models (Durech et al., 2015). 

After introducing the motivations and the different con¬ 
texts for calculating asteroid temperatures (§ 2), we pro¬ 
vide an overview of simple thermal models (§ 3) and of 
TPMs (§ 4). We describe data analysis techniques based 
on TPMs (§ 5), then we present the latest results and im¬ 
plications on the physics of asteroids (§ 6). In § 7, we dis¬ 
cuss temperature-induced surface changes on asteroids; see 
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also the chapter ’’Asteroid surface geophysics” by Murdoch 
et al. (2015). All used symbols are summarized in Tab. 1. 

Note that we do not discuss here the so-called asteroid 
thermal evolution models that are generally used to compute 
the temperature throughout the body as a function of time, 
typically taking into account internal heat sources such as 
the decay of the radiogenic Al. Such models allow one to 
estimate the degree of metamorphism, aqueous alteration, 
melting and differentiation that asteroids experienced dur¬ 
ing the early phases of the solar system formation (see Mc- 
Sween et al., 2002, for a review). 

2. MOTIVATIONS AND APPLICATIONS OF TPMs 

Thermophysical modeling of observations of asteroids 
in the thermal infrared (A > 4 /xm) is a powerful tech¬ 
nique to determine the values of physical parameters of as¬ 
teroids such as their sizes (e.g., Muller et al., 2014a), the 
thermal inertia and the roughness of their soils (e.g., Muller 
and Lagerros, 1998; Mueller, 2007; Delbo and Tanga, 2009; 
Matter et al., 2011; Rozitis and Green, 2014; Capria et al., 
2014) and in some particular cases also of their bulk density 
and their bulk porosity (Rozitis et al., 2013, 2014; Emery 
et al., 2014; Chesley et al., 2014) 

Knowledge of physical properties is crucial to under¬ 
stand asteroids: for instance, size information is fundamen¬ 
tal to constrain the asteroid size frequency distribution that 
informs us about the collisional evolution of these bodies 
(Bottke et al., 2005); is paramount for the study of aster¬ 
oid families, for the Earth-impact risk assessment of near- 
Earth asteroids (NEAs; see Harris et al., 2015, for a review), 
and for the development of asteroid space mission scenarios 
(§ 5.7). Accurate sizes are also a prerequisite to calculate 
the volumes of asteroids for which we know the mass, al¬ 
lowing us to derive the bulk density, which inform us about 
the internal structure of these bodies (e.g.. Carry, 2012). 

Thermal inertia, the resistance of a material to temper¬ 
ature change (§ 5.2), is a sensitive indicator for the prop¬ 
erties of the grainy soil (regolith, Murdoch et al., 2015) on 
asteroids, e.g., the typical grain size (Gundlach and Blum, 
2013) and their degree of cementation (Piqueux and Chris¬ 
tensen, 2009a,b) can be inferred from thermal-inertia mea¬ 
surements. In general, the regolith is what we study by 
means of remote-sensing observations. Understanding the 
regolith is therefore crucial to infer the nature of the under¬ 
lying body. Regolith informs us about the geological pro¬ 
cesses occurring on asteroids (Murdoch et al., 2015) such as 
impacts, micrometeoroid bombardment (Horz and Cintala, 
1997), and thermal cracking (Delbo et al., 2014). Regolith 
contains records of elements implanted by the solar wind 
and cosmic radiation, and therefore informs us about the 
sources of those materials (Lucey, 2006). Regolith porosity 
can shed light on the role of electrostatic and van-der-Waals 
forces acting on the surface of these bodies (e.g., Rozitis 
et al., 2014; Vernazza et al., 2012). 

Knowledge of surface temperatures is also essential for 
the design of the instruments and for the near-surface op¬ 


eration of space missions, as in the case of the sample- 
return missions Hayabusa-II and OSIRIS-REx of JAXA and 
NASA, respectively. In the future, knowledge of asteroid 
temperatures will be crucial for planning human interaction 
with asteroids. 

Another reason to model asteroid surface temperatures 
is that they affect its orbital and spin state evolution via 
the Yarkovsky and YORP effects, respectively (§ 5.8 and 
Vokrouhlicky et al., 2015). In particular, thermal inertia 
dictates the strength of the asteroid Yarkovsky effect. This 
influences the dispersion of members of asteroid families, 
the orbital evolution of potentially hazardous asteroids, and 
the delivery of 79 < 40 km asteroids and meteoroids from 
the main belt into dynamical resonance zones capable of 
transporting them to Earth-crossing orbits (see Vokrouh¬ 
licky et al., 2015, and the references therein). 

The YORP effect is believed to be shaping the distri¬ 
bution of rotation rates (Bottke et al., 2006) and spin vec¬ 
tor orientation (Vokrouhlicky et al., 2003; Hanus et al., 
2011, 2013); small gravitationally bound aggregates could 
be spun up so fast (Vokrouhlicky et al., 2015; Bottke et al., 
2006, and references therein) that they are forced to change 
shape and/or undergo mass shedding (Holsapple, 2010). 
Approximately 15% of near-Earth asteroids are observed to 
be binaries (Pravec et al., 2006), and YORP spin up is pro¬ 
posed as a viable formation mechanism (Walsh et al., 2008; 
Scheeres, 2007; Jacobson and Scheeres, 2011). 

A further motivation to apply TPM techniques is to con¬ 
strain the spin-axis orientation and the sense of rotation 
of asteroids (examples are 101955 Bennu and 2005 YUgs, 
Muller et al., 2012, 2013). Durech et al. (2015) describe 
how to use optical and thermal infrared data simultaneously 
to derive more reliable asteroid shapes and spin properties. 

The temperature and its evolution through the entire life 
of an asteroid can alter its surface composition and nature 
of the regolith (§ 7). Eor example, when the temperature 
rises above a certain threshold for a sustained period, certain 
volatiles can be lost via sublimation (Schorghofer, 2008; 
Capria et al., 2012), dehydration (Marchi et al., 2009), or 
desiccation (Delbo and Michel, 2011; Jewitt et al., 2015, 
and references therein). 

There can be pronounced and fast temperature variations 
between day and night. Modeling these temperature varia¬ 
tions is fundamental to studying the effect of thermal crack¬ 
ing of asteroid surface material (§ 7.1), which was found to 
be an important source of fresh regolith production (Delbo 
et al., 2014). 

3. SIMPLE THERMAL MODELS 

We start by introducing the near-Earth asteroid thermal 
model (NEATM, Harris, 1998) that is typically used where 
the data quality and/or the available knowledge about as¬ 
teroid shape and spin preclude the usage of TPMs. Typ¬ 
ically, the NEATM allows a robust estimation of asteroid 
diameter and albedo, but does not provide any direct infor¬ 
mation on thermal inertia or surface roughness (see Harris 
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and Lagerros, 2002, for a review). The recent large-scale 
thermal-emission surveys of asteroids and trans-Neptunian 
objects (see Mainzer et ah, 2015; Lellouch et ah, 2013, 
and references therein) typically use the NEATM in their 
data analysis, thereby establishing it as the de-facto default 
among the simple thermal models. The typical NEATM 
accuracy is 15% in diameter and roughly 30% in albedo 
(Harris, 2006). Other simple thermal models are the “Stan¬ 
dard Thermal Model” (STM; Lebofsky et ah, 1986), the 
“Isothermal Latitude Model” (ILM, also known as the East 
Rotating Model or ERM, Lebofsky and Spencer, 1989), and 
the night emission simulated thermal model (NESTM by 
Wolters and Green, 2009). The STM and the ILM, reviewed 
by, e.g., Harris and Lagerros (2002), have largely fallen out 
of use. 

The NEATM assumes that the asteroid has a spherical 
shape and does not directly account for thermal inertia nor 
surface roughness. The surface temperature is given by the 
instantaneous equilibrium with the insolation, which is pro¬ 
portional to the cosine of the angular distance between local 
zenith and the Sun and zero at the night side. The maximum 
temperature occurs at the subsolar point and it reads: 

(1-= ?7creT|s (1) 

(nomenclature is provided in Tab. 1). The parameter p was 
introduced in the STM of Lebofsky et al. (1986) as a means 
of changing the model temperature distribution to take ac¬ 
count of the observed enhancement of thermal emission at 
small solar phase angles due to surface roughness. This is 
known as the beaming effect. Eor this reason p is also called 
the beaming parameter. The 77 formalism, in the NEATM, 
allows a first-order description of the effect of a number 
of geometrical and physical parameters, in particular the 
thermal inertia and surface roughness on the spectral en¬ 
ergy distribution of an asteroid (Delbo et al., 2007). Eor 
a large thermal inertia, one would expect p-values signifi¬ 
cantly larger than unity (e.g., 1.5-3; with theoretical max¬ 
imum values around 3.5; Delbo et al., 2007), whereas for 
low thermal inertia p ~ 1 (for E = 0 and zero surface rough¬ 
ness). Roughness, on the other hand, tends to lower the 
value of rj (for observations at low or moderate phase an¬ 
gles). Eor instance, a value of p ~ 0.8 for a main belt as¬ 
teroid indicates that this body has low thermal inertia and 
significant roughness (with minimum theoretical values of 
0.6 - 0.7; Spencer, 1990; Delbo et al., 2007). We note, how¬ 
ever, that rj is not a physical property of an asteroid, as it can 
vary due to changing observing and illumination geometry, 
aspect angle, heliocentric distance of the body, phase angle 
and wavelength of observation. 

4. THERMOPHYSICAL MODELS 
4.1. Overview 


insulating surface that emitted thermal radiation in a non- 
Lambertian way (Pettit and Nicholson, 1930; Wesselink, 
1948a). Heat conduction and radiation scattering models 
of various rough surfaces were able to reproduce the lu¬ 
nar observations to a good degree (e.g.. Smith, 1967; Buhl 
et al., 1968b,a; Sexl et al., 1971; Winter and Krupp, 1971), 
and the derived thermal inertia and surface roughness val¬ 
ues matched in situ measurements by Apollo astronauts (see 
Rozitis and Green, 2011, and references therein). These 
early lunar models were adapted to general planetary bod¬ 
ies, albeit with an assumed spherical shape, by Spencer 
et al. (1989) and Spencer (1990). The most commonly 
used asteroid TPMs of Lagerros (1996a, 1997,1998), Delbo 
(2004), Mueller (2007), and Rozitis and Green (2011) are 
all based on Spencer et al. (1989) and Spencer (1990). Here 
we present the basic principles utilized in TPMs; for imple¬ 
mentation details, the reader is referred to the quoted works. 

All TPMs represent the global asteroid shape as a mesh 
of (triangular) facets (see Pig. 1) that rotates around a given 
spin vector with rotation period P. In general, utilized 
shape models are derived from radar observations, inversion 
of optical light-curves, in-situ spacecraft images, or stellar 
occultation timing (see Durech et al., 2015, for a review of 
asteroid shape modeling). If no shape model is available, 
one typically falls back to a sphere or an ellipsoid (e.g., 
Muller et al., 2013, 2014a; Emery et al., 2014). 

The goal is to calculate the thermal emission of each 
facet of the shape model at a given illumination and ob¬ 
servation geometry. To this end, the temperature of the sur¬ 
face and, in the presence of thermal inertia, the immediate 
subsurface need to be calculated. Generally, lateral heat 
conduction can be neglected as the shape model facets are 
much larger than the penetration depth of the diurnal heat 
wave (i.e., the thermal skin depth), and only ID heat con¬ 
duction perpendicular to and into the surface needs to be 
considered. Eor temperature T, time t, and depth z, ID 
heat conduction is described by: 


^dT d dT 


( 2 ) 


where k is the thermal conductivity, p is the material den¬ 
sity, and C is the heat capacity. If n is independent of depth 
(and, implicitly, temperature independent, see §5.3), Eq. 2 
reduces to the diffusion equation: 


dT K dT^ 

Ih ~ 


It is useful to define the thermal inertia E and the thermal 
skin depth Ig 

E = (4) 

Is = sJkPI2ttpC. (5) 


Different TPMs have been proposed to study the ther¬ 
mal emission of asteroids, comets, planets, and satellites. 
The first models were motivated by thermal observations 
of the lunar surface, which revealed an almost thermally 


These material properties are generally assumed to be 
constant with depth and temperature in asteroid TPMs, 
but varying properties has been considered in some Moon, 
Mars, planetary satellites, and asteroids models (e.g Giese 
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Table 1 

Nomenclature 


Symbol 

Quantity 

Unit 1 

1 Symbol 

Quantity 

Unit 

T 

Temperature 

K 

0 

Mean surface slope 

deg 

Tss 

Subsolar temperature 

K 

Is 

Thermal skin depth 

m 

a 

Stefan-Boltzmann constant 

(5.67051 10“®) Wm-2 K“^ 

P 

Rotation period 

s 

5© 

Solar constant at r=l au 

(1329) Wm“^ 

OJ 

= 2t7 IP 

s-i 

r 

Distance to the Sun 

au 

Xp 

ecliptic longitude of the pole 

deg 

r 

Vector to the Sun 

m 

/5p 

ecliptic latitude of the pole 

deg 

A 

Distance to the observer 

au 

<t>0 

Initial rotational phase at epoch 

deg 

€ 

Emissivity 

- 

a 

(Phase) angle between asteroid-sun-observer 

deg 

7] 

Beaming parameter 

- 

a 

Area of a facet 

vP 

K 

Thermal conductivity 

Wm“^ K“^ 

S 

Shadowing function 


c 

Heat capacity 

Jkg-^K”' 

P,I 

View factor 


r 

Thermal inertia 


Jvii) 

Visible radiosity 


0 

Thermal parameter 

- 

JiRiif) 

Infrared radiosity 


p 

Material density 

kgm-® 

h 

Local normal 

m 

H 

Absolute magnitude of the G system 


% 

Vector to the local facet 

m 

G 

Slope parameter of the H, G system 


3 

Vector to the remote facet 

m 

V 

Actual magntude in the V-band 


1C 

Crater opening angle 

deg 

D 

Diameter 

m (or km) 

PC or/ 

Area density of craters 


A 

Bolometric Bond albedo 

- 


Emission angle 

rad 

Pv 

Geometric visible albedo 

- 

/a(t) 

Infrared flux 


z 

Depth in the subsoil 

m 

A 

Wavelength 

/j.m 

t 

Time 

s 




rp 

Pore radius of regolith 

m 





and Kuehrt, 1990; Urquhart and Jakosky, 1997; Piqueux 
and Christensen, 2011; Keihm, 1984; Keihm et ak, 2012; 
Capria et ak, 2014, see also section 5.3). 

TPM implementations typically employ dimensionless 
time and depth variables; t = 2Trt/P and Z = z/ls- Then, 
the only remaining free parameter is the dimensionless ther¬ 
mal parameter (0 = ry/uj/eaTgg, Spencer, 1990) describ¬ 
ing the combined effect of thermal inertia, rotation period, 
and heat emission into space on the surface temperature dis¬ 
tribution (see Fig. 2). 

A numerical finite-difference technique is used to solve 
the ID heat conduction equation, and an iterative technique 
is used to solve the surface boundary condition. This re¬ 
quires a suitable number of time and depth steps to fully 
resolve the temperature variations and to ensure model sta¬ 
bility (typically, at least 360 time steps and 40 depth steps 
over 8 thermal skin depths are required). TPMs are run 
until specified convergence criteria are met (e.g., until tem¬ 
perature variations between successive model iterations are 
below a specified level) and/or until a specified number of 
model iterations have been made. 

For applications such as the study of the sublimation of 
water ice from the shallow subsurface of asteroids (e.g., 
the Main Belt Comets or 24 Themis) the heat conduction 
equation must be coupled with a gas diffusion equation 
(Schorghofer, 2008; Capria et ak, 2012; Prialnik and Rosen¬ 
berg, 2009). See also Huebner et ak (2006) for a review. 

The ID heat conduction equation is solved with inter¬ 
nal and surface boundary conditions to ensure conservation 
of energy. Since the amplitude of subsurface temperature 
variations decreases exponentially with depth, an internal 
boundary condition is required to give zero temperature gra¬ 


dient at a specified large depth 


Ih J 


= 0 . 




( 6 ) 


A typical surface boundary condition for a facet at point i 
with respect to the asteroid origin, at point r with respect to 
the Sun, and with surface normal h is then given by 


eaT‘^{i, t) — 


dT{i, t) 
dz 


z=0 


J Jv{j)Fjjda'+ 

ecr(l-e) J JmijjFrjda' 


(7) 


The left-hand side of Eq. (7) gives the thermal energy ra¬ 
diated to space and the heat conducted into the subsurface, 
and the right-hand side gives the input radiation from three 
different sources; direct solar radiation, multiply scattered 
solar radiation (i.e., self-illumination), and reabsorbed ther¬ 
mal radiation (i.e., self-irradiation). The two last compo¬ 
nents are also known as mutual heating (see Fig. 3). 

The amount of solar radiation absorbed by a facet de¬ 
pends on the Bond albedo A and any shadows projected on 
it, which is dictated by (i.e., S{f,i) = 1 or 0, de¬ 

pending on the presence or absence of a shadow). Projected 
shadows occur on globally non-convex shapes only, which 
can be determined by ray-triangle intersection tests of the 
solar illumination ray (e.g., Rozitis and Green, 2011) or by 
local horizon mapping (e.g., Statler, 2009). Related to shad¬ 
owing are the self-heating effects arising from interfacing 
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Fig. 1. — (a) example of a triangulated 3D shape model 
as typically used in TPMs (asteroid (2063) Bacchus from 
http://echo.jpl.nasa.gov/asteroids/shapes/shapes.html). Tempera¬ 
tures are color coded: white corresponds to the maximum and 
dark-gray corresponds to minimum temperature. Three different 
roughness models are sketched in the bottom of the figure: (b) 
hemispherical section craters; (c) Gaussian surface; (d) fractal sur¬ 
face. Sub-figures b and d are adapted from Davidsson et al. (2015), 
c is from Rozitis and Green (2011). 

facets, which tend to reduce the temperature contrast pro¬ 
duced inside concavities. The problem here is to determine 
which facets see other facets, and to calculate the amount 
of radiation exchanged between them. The former can be 
determined by ray-triangle intersection tests again, and the 
latter can be solved using view factors. The view factor 
Fr y is defined as the fraction of the radiative energy leaving 
the local facet i that is received by the remote facet J as¬ 
suming Lambertian emission (Lagerros, 1998). and 

JiR^j) are then the visible and thermal-infrared radiosities 
of remote facet j. Either single or multiple scattering can be 
taken into account, and the latter can be efficiently solved 
using Gauss-Seidel iterations (Vasavada et al., 1999). Most 
TPMs neglect shadowing and self-heating effects resulting 
from the global shape for simplicity, but they can be sig¬ 
nificant on asteroids with large shape concavities (e.g., the 
south pole of (6489) Golevka: Rozitis and Green, 2013). 

4.2. Modeling asteroid thermal emission 

Once the surface temperature distribution across an as¬ 
teroid surface has been computed, the emission spectrum 
(Fig. 4) at a given observation geometry and a specified time 
can be calculated. The monochromatic flux density can also 
be calculated at wavelengths of interest. When these model 
fluxes are plotted as a function of the asteroid rotational 



Time (h) 

Fig. 2.— Synthetic diurnal temperature curves on the equator of 
a model asteroid for different values of thermal inertia (in units 
of J m^^s '^K ’). Increasing thermal inertia smooths temperature 
contrasts and causes the temperature peak to occur after the insola¬ 
tion peak at 3 h. The asteroid is situated at a heliocentric distance 
of r = 1.1 au, has a spin period of 6 h, a Bond albedo of A = 0.1, 
and its spin axis is perpendicular to the orbital plane. 

phase, one obtains the so-called thermal lightcurves (e.g. 
Fig. 5), which can be used to test the fidelity of shape and 
albedo models typically used as input in the TPM. 

When the temperature for a facet is known, the inten¬ 
sity I\ (t) at which it emits at wavelength A is given by the 
Planck function. Assuming Lambertian emission, the spec¬ 
tral flux of the facet seen by an observer is then 

fxir) = cos cj) (8) 

where a is area of the facet, A is the distance to the ob¬ 
server, and (j) is the emission angle. The total observed 
flux is obtained by summing the thermal fluxes of all vis¬ 
ible shape model facets including any contributions from 
surface roughness elements contained within them. For 
disk-integrated measurements, this summation is performed 
across the entire visible side of the asteroid, whilst for spa¬ 
tially resolved measurements it is summed across facets 
contained within the detector pixel’s held of view. 

The assumption of Lambertian emission depends on no 
directionality induced by surface irregularities at scales 
below the thermal skin depth. Davidsson and Rickman 
(2014) show that surface roughness at sub-thermal-skin- 
depth scales is quasi-isothermal and is therefore not likely 
to deviate from Lambertian emission overall. However, 
radiative transfer processes between the regolith grains 
could contribute up to 20% of the observed beaming effects 
(Hapke, 1996). Rozitis and Green (2011, 2012) investi¬ 
gated combined microscopic (regolith grain induced) and 
macroscopic (surface roughness induced) beaming effects, 
and demonstrated that the macroscopic effects dominated 
overall. This was previously found to be the case in di¬ 
rectional thermal emission measurements of lava flows on 
Earth (Jakosky et al., 1990). 
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Fig. 3.— Diagram illustrating the energy balance and radiation 
transfer between facets (copied from Rozitis and Green, 2011). 
The terms Fsun, Fscat, Frad, k{dT/dx) and eaT* are the di¬ 
rect sunlight, multiply scattered sunlight, reabsorbed thermal radi¬ 
ation, conducted heat and thermal radiation lost to space, respec¬ 
tively. 

As wavelengths increase to the submillimeter range and 
above, asteroid regolith becomes increasingly transparent 
and the observed flux is integrated over increasing depths 
(Chamberlain et al., 2009; Keihm et al., 2013). Model¬ 
ing such fluxes with typical thermal models (which derive 
fluxes from surface temperatures, only) requires a signifi¬ 
cant reduction in effective spectral emissivity. For example, 
3.2 mm flux measurements of (4) Vesta require an emissiv¬ 
ity of ~0.6 to match model predictions (Muller and Barnes, 
2007). The reduction in emissivity can be explained by 
lower subsurface emission temperatures (Lagerros, 1996b) 
and by different subsurface scattering processes dependent 
on grain size (Redman et al., 1992; Muller and Lagerros, 
1998). Keihm et al. (2012,2013) attribute the reduced emis¬ 
sivity at submm/mm wavelengths to a higher thermal inertia 
value of the subsurface layers. Reduction in emissivity has 
also been determined at wavelengths shorter than 4.9 /im 
for disk-resolved regions of (4) Vesta (Tosi et al., 2014). 

4.3. Surface Roughness 

Roughness causes an asteroid surface to thermally emit 
in a non-Lambertian way with a tendency to reradiate the 
absorbed solar radiation back towards the Sun, an effect 
known as thermal infrared beaming (Lagerros, 1998; Rozi¬ 
tis and Green, 2011). It is thought to be the result of two 
different processes: a rough surface will have elements ori¬ 
entated towards the Sun that become significantly hotter 
than a flat surface, and multiple scattering of radiation be¬ 
tween rough surface elements increases the total amount of 
solar radiation absorbed by the surface. The relevant size 
scale ranges from the thermal skin depth to the linear size 
of the facets in the shape model. It is included in thermo¬ 
physical models by typically modeling an areal fractional 



r.I.i... 

0 10 20 30 40 

Wavelength (micron) 


Fig. 4.— (a) Example of SED calculated from a TPM and 
the NEATM compared to a Spitzer spectrum of (87) Sylvia, (b) 
Spectral emissivity derived from the above: data divided by the 
NEATM continuum. Eigure from Marchis et al. (2012). 

coverage (/) of spherical-section craters (of opening an¬ 
gle 7 c) within each shape model facet. Other more com¬ 
plex forms have been considered, such as Gaussian random 
(Lagerros, 1998) or fractal (Groussin et al., 2013) surfaces 
or parallel sinusoidal trenches (see sketch of Fig. 1), but the 
spherical-section crater produces similar results (in terms of 
the disk-integrated beaming effect Lagerros, 1998) and ac¬ 
curately reproduces the directionality of the lunar thermal 
infrared beaming effect (Rozitis and Green, 2011). How¬ 
ever, it has been shown that the thermal emission depends 
also on roughness type in addition to roughness level, for 
disk resolved data Davidsson et al. (2015). 

Spherical-section craters are typically implemented, as 
the required shadowing and view-factor calculations can 
be performed analytically (Emery et al., 1998; Lagerros, 
1998). Heat conduction can be included by dividing the 
crater into several tens of surface elements where the same 
equations listed above can be applied. Alternatively, the 
temperature distribution within the crater resulting from 
heat conduction, rcjater(r), can be approximated using 

Though (F) _ rsmooth(r) 

Trough ( 0 ) Tsmooth ( 0 ) 

where Trough (0) can be calculated analytically assuming in¬ 
stantaneous equilibrium (Lagerros, 1998). rsniooth(O) and 
Tsmooth(r) are the corresponding smooth-surface tempera¬ 
tures, which can be calculated exactly. This approximation 
is computationally much cheaper than the full implemen¬ 
tation. However, it does not work on the night side of the 
asteroid and temperature ratios diverge near the terminator 
(Mueller, 2007). An even simpler alternative is to multi¬ 
ply the smooth-surface temperatures by a NEATM-like r] 
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Fig. 5.— Example of a TPM generated thermal lightcurve 
(dashed line) and real data for (101955) Bennu. From Emery et al. 
(2014). 

value (e.g. Groussin et al., 2011). Whilst this alternative 
might produce the correct disk-integrated color tempera¬ 
ture of the asteroid, it does not reproduce the directionality 
of the beaming effect. Indeed, roughness models predict a 
limb-brightening effect (Rozitis and Green, 2011), which is 
seen in spatially-resolved measurements of (21) Lutetia by 
Rosetta (Keihm et al., 2012). 

The above implementations neglect lateral heat conduc¬ 
tion, although the spatial scales representing surface rough¬ 
ness can, in some cases, become comparable to the thermal 
skin depth. Modeling of 3D heat conduction inside rocks 
the size of the thermal skin depth has demonstrated that 
their western sides (for a prograde rotator; eastern sides for 
a retrograde rotator) are generally warmer than their east¬ 
ern sides, which could result in a tangential-YORP effect 
that predominantly spins up asteroids (Golubov and Krugly, 
2012). Other than this, it appears that the ID heat conduc¬ 
tion approximation still produces satisfactory results. 

In thermophysical models, the degree of surface rough¬ 
ness can be quantified in terms of the Hapke mean surface 
slope 

2 

tan0=— / a{6)taYi9 dO (10) 

TT Jo 

where 6 is the angle of a given facet from horizontal, and 
a{6) is the distribution of surface slopes (Hapke, 1984). Al¬ 
ternatively, it can be measured in terms of the RMS sur¬ 
face slope (Spencer, 1990). This then allows comparison of 
results derived using different surface roughness represen¬ 
tations (e.g., craters of different opening angles and frac¬ 
tional coverages, or different Gaussian random surfaces: 
Davidsson et al., 2015), and comparison against rough¬ 
ness measured by other means. It has been demonstrated 
that different roughness representations produce similar de¬ 
grees of thermal infrared beaming when they have the same 
degree of roughness measured in terms of these values 
(Spencer, 1990; Emery et al., 1998; Lagerros, 1998; Rozitis 


and Green, 2011). 

5. DATA ANALYSIS USING A TPM 
5.1. Thermal infrared spectro-photometry 

Physical properties than can be derived from TPM fits 
to disk-integrated thermal observations include the diam¬ 
eter, geometric albedo, thermal inertia and roughness. In 
practically all cases, the absolute visual magnitude H is 
known, establishing a link between D and pv and reduc¬ 
ing the number of TPM fit parameters by one: 

D(km) = 1329(11) 

(Fowler and Chillemi, 1992; Vilenius et al., 2012). Fre¬ 
quently, the rotational phase during the thermal observa¬ 
tions is not sufficiently well known and has to be fitted to 
the thermal data (e.g., Harris et al., 2005; All-Lagoa et al., 
2014). In some cases, TPMs can be used to constrain the 
orientation of the spin vector of an asteroid, with Ap and jjp 
treated as free parameters (as demonstrated e.g., by Muller 
et al., 2013, 2012, note that in the case of 101955 Bennu the 
radar-constrained pole solution was not yet known). More¬ 
over, Muller et al. (2014b) successfully performed a TPM 
analysis of an asteroid (99942 Apophis) in a non-principal 
axis rotation state for the first time. 

The thermal effects of thermal inertia and surface rough¬ 
ness are difficult to tell apart. A commonly used approach is 
to use four different roughness models corresponding to no, 
low, medium, and high roughness, with each model leading 
to a different thermal-inertia fit (Mueller, 2007; Delbo and 
Tanga, 2009); frequently, the scatter between these four so¬ 
lutions accounts for the bulk of the uncertainty in thermal 
inertia. However, in some lucky cases, data do allow the 
effects of roughness and thermal inertia to be disentangled. 
This requires good wavelength coverage straddling the ther¬ 
mal emission peak and good coverage in solar phase angle, 
such that both the morning and afternoon sides of the aster¬ 
oid are seen. See Fig. 6 for an illustration. 

The best-fitting model parameters are those that mini¬ 
mize Their uncertainty range is spanned by the values 
that lead to within a specified threshold of the best fit, 
depending on the number of free fit parameters. Ideally, the 
reduced x^ of the best fit should be around unity. How¬ 
ever, due to systematic uncertainties introduced in thermal 
infrared observations (e.g., flux offsets between different in¬ 
struments), and/or in the thermophysical modeling, it is not 
uncommon to get large reduced x^ values. Large x^ val¬ 
ues are also obtained when the assumed shape model dif¬ 
fers significantly from the asteroid’s true shape (Rozitis and 
Green, 2014). In particular, if the spatial extent of the shape 
model’s z-axis is wrong, this can lead to diameter determi¬ 
nations that are inconsistent with radar observations (e.g., 
for 2002 NY40 and (308635) 2005 YU55 in Muller et al., 
2004, 2013, respectively), and/or two different thermal in¬ 
ertia determinations (e.g., the two different results produced 
for (101955) Bennu by Emery et al., 2014; Muller et al.. 
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Fig. 6.— (21) Lutetia: a TPM fit that allows surface rough¬ 
ness to be constrained. The quantity p here is not the bulk 
density of the body, but it is the r.m.s. of the slopes on 
the surface. It is related to the ratio between the diameter 
and the depth of spherical section craters (Lagerros, 1998) 
in this particular case. / is the areal fraction of each facet 
covered with craters. From O’Rourke et al. (2012). 

2012). In some works, the asteroid shape model has also 
been optimized during the thermophysical fitting to resolve 
inconsistencies with radar observations (e.g., (1862) Apollo 
and (1620) Geographos in Rozitis et al., 2013; Rozitis and 
Green, 2014, respectively). 

We remind the reader that the accuracy of the physical 
properties (in particular the value of F) of asteroids derived 
from TPM depends on the quality of the thermal infrared 
data, coverage in wavelength, phase, rotational, and aspect 
angle. The accuracy of the shape model and of the H and 
G values are also important (see e.g. Rozitis and Green, 
2014). The derived thermal inertia value often depends on 
the assumed degree of roughness and it is usually affected 
by large errors (e.g. 50 or 100 %, see Tab. 2). Care must be 
used in accepting TPM solutions purely based on the good¬ 
ness of fit (e.g. the value of the x^), as they can be domi¬ 
nated by one or few measurements with unreliable small er¬ 
rors or calibration offsets between measurements from dif¬ 
ferent sources. 

5.2. Thermal Inertia and Thermal Conductivity 

As asteroids rotate, the day-night cycle causes cyclic 
temperature variations that are controlled by the thermal in¬ 
ertia (defined by Eq. 4) of the soil and the rotation rate of the 
body. In the limit of vanishing thermal inertia, the surface 
temperature would be in instantaneous equilibrium with the 
incoming flux, depending only on the solar incidence an¬ 
gle (as long as self heating can be neglected); surface tem¬ 
peratures would peak at local noon and would be zero at 
night. In reality, thermal conduction into and from the sub¬ 
soil causes a certain thermal memory, referred to as ther¬ 
mal inertia. This smoothens the diurnal temperature profile. 


leads to non-zero night-side temperatures, and causes the 
surface temperature to peak on the afternoon side, as shown 
in Fig. 2, thereby causing the Yarkovsky effect (§ 5.8). 

The mass density p, the specific heat capacity c, and the 
thermal conductivity k, and correspondingly F itself, must 
be thought of as effective values representative of the depth 
range sampled by the heat wave, which is typically in the 
few-cm range. In turn, thermal-inertia values inform us 
about the physical properties of the top few centimeters of 
the surface, not about bulk properties of the object. 

As will be discussed below (§ 6.5), p and c of an asteroid 
surface can plausibly vary within a factor of several, while 
plausible values of k span a range of more than 4 orders of 
magnitude. It is therefore not unjustified to convert from 
r to K and back using reference values for p and c (note 
that Yarkovsky/YORP models tend to phrase the thermal- 
conduction problem in units of n, while TPMs tend to be 
formulated in units of F, which is the observable quantity). 

Importantly, the k of finely powdered lunar regolith is 3 
orders of magnitude lower than that of compact rock (com¬ 
pact metal is even more conductive by another order of 
magnitude). This is because radiative thermal conduction 
between regolith grains is significantly less efficient than 
phononic heat transfer within a grain. A fine regolith, an 
aggregate of very small grains, is a poor thermal conduc¬ 
tor and displays a low F. Thermal inertia can therefore 
be used to infer the presence or absence of thermally in¬ 
sulating powdered surface material. In extension, thermal 
inertia can be used as a proxy of regolith grain size. The 
required calibration under Mars conditions (where the ten¬ 
uous atmosphere enhances thermal conduction within pores 
compared to a vacuum) was obtained by Presley and Chris¬ 
tensen (1997a,b) and used in the analysis of thermal-inertia 
maps of Mars (see Mellon et al., 2000; Putzig and Mel¬ 
lon, 2007). Similar progress in asteroid science was slowed 
down by the lack of corresponding laboratory measure¬ 
ments under vacuum conditions (but see below for recent 
lab measurements of meteorites). However, Gundlach and 
Blum (2013) provided a calibration relation based on heat- 
transfer modeling in a granular medium. 

5.3. Temperature dependence of thermal inertia 

Thermal inertia is a function of temperature (Keihm, 
1984), chiefly because the thermal conductivity is. In gen¬ 
eral, for a lunar-like regolith the thermal conductivity is 
given by: 

«: = Kf, + AaVpT^ (12) 

where Kb is the solid-state thermal conductivity (heat con¬ 
duction by phonons) and Vp is the radius of the pores of the 
regolith. The term proportional to is due to the heat con¬ 
duction by photons. Equation 12 is often written in the form 
K = Ko(l + xT^) (6-g-> Vasavada et al., 1999). Note that 
Kb is itself a function of T (Opeil et al., 2010). There is ex¬ 
tensive literature on the T-dependence of the conductivity 
of lunar regolith (e.g., Vasavada et al., 1999, and references 
therein). A theoretical description of the temperature de- 







pendence of k in regoliths is given by Gundlach and Blum 
(2013). 

Asteroid TPMs typically neglect the temperature depen¬ 
dence of r. This is uncritical for typical remote observa¬ 
tions, which are dominated by the warm sunlit hemisphere 
(see Fig. 7 and Capria et al., 2014; Vasavada et al., 2012). In 
the analysis of highly spatially resolved observations, how¬ 
ever, the temperature dependence must be considered, cer¬ 
tainly when analyzing night-time observations on low-F as¬ 
teroids. Note that for temperature-dependent k, Eq. 2 must 
be used instead of Eq. 3 (see, e.g., Capria et al., 2014). 



Eig. 7.— Diurnal temperature curves at the equator of an asteroid 
with A=0.1, P=6h, e=1.0 at a heliocentric distance of 1.0 au. Solid 
curve: constant thermal conductivity h;=0.02 W K * m '. Dashed- 
curved: temperature-dependent heat conductivity k = 10“^(1 -f 
0.5(T/250)^) W K ‘ m '. 

Caution must be exercised when comparing thermal- 
inertia results obtained at different heliocentric distances r, 
i.e., at different temperatures. All other things being equal, 
oc r“^. Assuming that the term dominates in Eq. 
12, the thermal inertia of a test object scales with (see also 
Mueller et al., 2010): 

r oc Vk oc T3/2 oc r-3/^. (13) 

5.4. Binary Asteroid TPM 

A rather direct determination of thermal inertia can be 
obtained by observing the thermal response to eclipses and 
their aftermath, allowing one to see temperature changes in 
real time. Such observations have been carried out for plan¬ 
etary satellites such as the Galilean satellites (Morrison and 
Cruikshank, 1973), and our Moon (Pettit, 1940; Shorthill, 
1973; Lawson et al., 2003; Lucey, 2000; Eountain et al., 
1976). Mueller et al. (2010) report the first thermal obser¬ 
vations of eclipse events in a binary Trojan asteroid system, 
(617) Patroclus, where one component casts shadow on the 
other while not blocking the line of sight toward the ob¬ 
server. 

The thermophysical modeling of eclipse events is rela¬ 
tively straightforward, assuming the system is in a tidally 
locked rotation state typical of evolved binary systems. In 


that case, the components’ spin rates match one another, 
and their spin axes are aligned with that of the mutual or¬ 
bit. The system is therefore at rest in a co-rotating frame 
and can be modeled like a single object with a non-convex 
(disjoint!) global shape. Eclipse effects are fully captured, 
provided that shadowing between facets is accounted for. 
The two hemispheres that face one another can, in princi¬ 
ple, exchange heat radiatively. This is negligible for typical 
binary systems, however. 

As discussed above, the thermal effects of roughness and 
thermal inertia can be hard to disentangle. In the case of 
eclipse measurements, which happen at essentially constant 
solar phase angle, the effect of surface roughness is much 
less of a confounding factor. This is because the variation in 
the thermal signal is dominated by the temperature change 
induced by the passing shadow, which is a strong function 
of thermal inertia. 

It must be kept in mind that the duration of an eclipse 
event is short compared to the rotation period. The eclipse- 
induced heat wave therefore probes the subsoil less deeply 
than the diurnal heat wave does (the typical heat penetration 
depth is given by Eq. 5 with P equal to the duration of the 
eclipse event). A depth dependence of thermal inertia (see 
§ 6.6) could manifest itself in different thermal-inertia de¬ 
terminations using the two different measurement methods. 

5.5. Thermal-infrared interferometry 

Interferometric observations of asteroids in the ther¬ 
mal infrared measure the spatial distribution of the ther¬ 
mal emission along different directions on the plane of sky, 
thereby constraining the distribution in surface temperature 
and hence thermal inertia and roughness. Provided the as¬ 
teroid shape is known, interferometry can be used to break 
the aforementioned degeneracy between thermal inertia and 
roughness from a single-epoch observation (Matter et al., 
2011, 2013). Interferometry also allows a precise determi¬ 
nation of the size of an asteroid (Delbo et al., 2009). 

Spatial resolutions between 20 and 200 milli-arcseconds 
can be obtained from the ground (see Durech et al., 2015, 
for a review and future perspectives of the application of 
this technique). 

While for the determination of asteroid sizes and shapes 
from interferometric observations in the thermal infrared, 
simple thermal models can be used (Delbo et al., 2009; 
Carry et al., 2015), a TPM was utilized to calculate interfer¬ 
ometric visibilities of asteroids in the thermal infrared for 
the observations of (41) Daphne (Matter et al., 2011) and 
(16) Psyche (Matter et al., 2013). 

Mid-infrared interferometric instruments measure the to¬ 
tal flux and the visibility of a source, the latter being re¬ 
lated to the intensity of the Eourier Transform (ET) of the 
spatial flux distribution along the interferometer’s baseline 
projected on the plane of sky. Thus, the data analysis proce¬ 
dure consisted in generating images of the thermal infrared 
emission of the asteroids at different wavelengths as viewed 
by the interferometer and then in obtaining the model visi- 


9 






bility and flux for each image. The former is related to the 
FT of the image, the latter is simply the sum of the pix¬ 
els. The free parameters of the TPM (size, thermal inertia 
value and roughness) are adjusted in order to minimize the 
distance between the disk integrated flux and visibility of 
the model, and the corresponding observed quantities (see 
Matter et al., 2011, 2013; Durech et al., 2015, for further 
information). Some results from these observational pro¬ 
grams are discussed in S 6 and in the chapter by Durech 
et al. (2015). 

5.6. Disk-resolved data and retrieval of temperatures 

The availability of disk-resolved thermal-infrared obser¬ 
vation has been increasingly steadily over the years: the 
ESA mission Rosetta performed two successful flybys to 
the asteroids (2867) Steins (in 2008) and (21) Lutetia (in 
2010) (Barucci et al., 2015). In 2011, the NASA mission 
Dawn began its one-year orbiting of (4) Vesta (Russell et al., 
2012); Dawn is reaching (1) Ceres at the time of writing. 
JAXA’s Hayabusa II sample-return mission will map its tar¬ 
get asteroid (162173) 1999 JU 3 using the thermal infrared 
imager (TIR) aboard the spacecraft (Okada et al., 2013); 
NASA’s OSIRIS-REx mission and its thermal spectrome¬ 
ter OTES will do likewise for its target asteroid (101955) 
Bennu (in 2018-2019). These data can be used to derive 
surface-temperature maps, from which maps of thermal in¬ 
ertia and roughness can be derived. 

Three different methods are used to measure surface 
temperatures from orbiting spacecraft: bolometry, mid- 
infrared spectroscopy, and near-infrared spectroscopy. In 
the following, we will elaborate on the challenges posed by 
these different methods, and on their dependence on spec¬ 
tral features, surface roughness, illumination geometry, and 
viewing geometry. 

Bolometers measure thermal flux within a broad band¬ 
pass in the infrared, approximating the integral of the 
Planck function, U = aT^ (e-g-, Kieffer et al., 1977; Paige 
et al., 2010). The temperature derived in this way (effective 
temperature) is directly relevant to the energy balance 
on the surface. Since the bolometric flux is spectrally 
integrated, the resulting temperature is fairly insensitive 
to spectral emissivity variations, as long as the bolomet¬ 
ric emissivity (weighted spectrally averaged emissivity) is 
known or can be reasonably approximated. 

Temperatures derived from mid-infrared spectrometry, 
on the other hand, are typically brightness temperatures, 
i.e., the temperature of a black body emitting at the wave¬ 
length in question. It is generally assumed that at some 
wavelength, the spectral emissivity is very close to 1 . 0 , and 
the brightness temperature at this wavelength is taken as the 
surface temperature. 

Spacecraft sent to asteroids (and/or comets) have more 
commonly been instrumented with near-infrared spectrom¬ 
eters (e.g., A < 5 iim) rather than mid-infrared spectrome¬ 
ters. The long-wavelength ends of these spectrometers of¬ 
ten extend into the range where thermal emission dominates 


the measured flux (for the daytime surface temperatures of 
most asteroids). At these wavelengths, one cannot assume 
that the emissivity is close to 1.0. It is therefore not prac¬ 
tical to derive brightness temperatures. Instead, the color 
temperature is derived, that is the temperature of a black 
body that emits with the same spectral shape. Such deriva¬ 
tions have to separate temperature from spectral emissivity. 
The problem is under-constrained (Nh- 1 unknowns, but only 
N data points), so there is no deterministic solution. Spec¬ 
tral emissivities for fine-grained silicates trend in the same 
direction as the blackbody curve, so it would be very easy 
to mistake spectral emissivity variations for different tem¬ 
peratures. The most statistically rigorous approach that has 
been applied to separating temperature and spectral emis¬ 
sivity in the 3-5 /rm region is that of Keihm et al. (2012) 
and Tosi et al. (2014) for Rosetta/VIRTIS data of Lutetia 
and Dawn/VIR data of Vesta. 

Temperatures thus measured represent an average tem¬ 
perature in the field of view of a given pixel: illuminated hot 
zones and shadowed colder parts will both contribute. They 
do not directly correspond to a physical temperature of the 
soil; rather, they depend sensitively on the observation and 
illumination geometry (see Rozitis and Green, 2011, in par¬ 
ticular their Eig. 9), especially in the case of large illumina¬ 
tion angles. 

Microwave spectrometers such as MIRO (Gulkis et al., 
2007) can provide both day and nightside thermal flux mea¬ 
surements. At sub-mm, mm, and longer wavelengths, aster¬ 
oid soils become moderately transparent. Subsurface lay¬ 
ers contribute significantly to the observable thermal emis¬ 
sion, thus providing information on the subsurface temper¬ 
ature. Observable fluxes depend on the subsurface temper¬ 
ature profile, weighted by the wavelength-dependent elec¬ 
trical skin depths, so both a thermal and an electrical model 
are required to interpret such data (Keihm et al., 2012). 

We remind here that thermal infrared fluxes should be 
used as input data for TPMs and not (effective, color, or 
brightness) temperatures derived from radiometric meth¬ 
ods, because of their dependence on illumination and ob¬ 
servation angles ! 

5.7. Sample-return missions 

Space agencies across the planet are developing space 
missions to asteroids, notably sample-return missions to 
primitive (C and B type) near-Earth asteroids: Hayabusa- 
2, was launched by JAXA towards (162173) 1999 JU 3 on 
December 3, 2014, and OSIRIS-REx is to be launched by 
NASA in 2016 (Lauretta et al., 2012). A good understand¬ 
ing of the expected thermal environment, which is governed 
by thermal inertia, is a key factor in planning spacecraft op¬ 
erations on or near asteroid surfaces. E.g., OSIRIS-REx is 
constrained to sampling a regolith not hotter than 350 K, 
severely constraining the choice of the latitude of the sam¬ 
ple selection area on the body, the local time, and the arrival 
date on the asteroid. 

Both Hayabusa-2 and OSIRIS-REx are required to take 
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regolith samples from the asteroid surface back to Earth. 
Obviously, this requires that regolith be present in the first 
place, which needs to be ascertained by means of ground- 
based thermal-inertia measurements. The sampling mecha¬ 
nism of OSIRIS-REx, in particular, requires relatively fine 
(cm-sized or smaller) regolith. 

5.8. Accurate Yarkovsky and YORP modeling from 
TPMs 

Scattered and thermally emitted photons carry momen¬ 
tum. Any asymmetry in the distribution of outgoing pho¬ 
tons can, after averaging over an orbital period, impart a net 
recoil force (Yarkovsky effect) and/or a net torque (YORP 
effect) on the asteroid. Both effects are more noticeable as 
the object gets smaller. Eor small enough objects, the orbits 
can be significantly affected by the Yarkovsky effect, and 
their rotation state by YORP (Bottke et ak, 2006; Vokrouh- 
licky et ak, 2015). 

The strength of the Yarkovsky effect is strongly influ¬ 
enced by thermal inertia (Bottke et ak, 2006, and references 
therein) and by the degree of surface roughness (Rozitis 
and Green, 2012). However, the strength and sign of the 
YORP rotational acceleration on an asteroid is independent 
of thermal inertia (Capek and Vokrouhlicky, 2004), but it 
is highly sensitive to the shadowing (Breiter et ak, 2009), 
self-heating (Rozitis and Green, 2013), and surface rough¬ 
ness effects (Rozitis and Green, 2013) that are incorporated 
in thermophysical models. 

Accurate calculations of the instantaneous recoil forces 
and torques require an accurate calculation of surface tem¬ 
peratures as afforded by TPMs; Rozitis and Green (2012, 
2013) report on such models. Other than on thermal in¬ 
ertia, the Yarkovsky-induced orbital drift depends on the 
bulk mass density. Therefore, Yarkovsky measurements 
combined with thermal-inertia measurements can be used 
to infer the elusive mass density (Mommert et ak, 2014b,a; 
Rozitis and Green, 2014; Rozitis et ak, 2014, 2013). In the 
case of (101955) Bennu, the uncertainties in published val¬ 
ues of thermal inertia (Emery et ak, 2014) and measured 
Yarkovsky drift (Chesley et ak, 2014) are so small that the 
accuracy of the inferred mass density rivals that of the ex¬ 
pected in-situ spacecraft result (1260 ± 70 kg m“^, i.e., a 
nominal uncertainty of only 6%). Rozitis et ak (2014) de¬ 
rived the bulk density of (29075) 1950 DA and used it to 
reveal the presence of cohesive forces stabilizing the object 
against the centrifugal force. 

In turn, the measured Yarkovsky drift can be used to infer 
constraints on thermal inertia. This was first done by Ches¬ 
ley et ak (2003) studying the Yarkovsky effect on (6489) 
Golevka and by Bottke et ak (2001) studying the Koro- 
nis family in the main asteroid belt; both studies revealed 
thermal inertias consistent with expectations based on the 
observed correlation between thermal inertia and diameter 
(see § 6.2). 

Whilst the YORP effect is highly sensitive to small-scale 
uncertainties in an asteroid’s shape model (Statler, 2009), it 


can be used to place constraints on the internal bulk density 
distribution of an asteroid (Scheeres and Gaskell, 2008). 
Eor instance, Lowry et ak (2014) explain the YORP de¬ 
tection on (25143) Itokawa, which was opposite in sign to 
that predicted, by Itokawa’s two lobes having substantially 
different bulk densities. However, unaccounted lateral heat 
conduction in thermal skin depth sized rocks could also ex¬ 
plain, at least partially, this opposite sign result (Golubov 
and Krugly, 2012). 

6. LATEST RESULTS FROM TPMs 

In the Asteroid III era, thermal properties were known 
for only a few asteroids, i.e., (1) Ceres, (2) Pallas, (3) 
Juno, (4) Vesta, (532) Herculina from Muller and Lagerros 
(2002), and (433) Eros from Lebofsky and Rieke (1979). 
Since then, the number of asteroids with known thermal 
properties has increased steadily. We count 59 minor bodies 
with known value of E (see Tab. 2). Of these, 16 are near- 
Earth asteroids (NEAs), 27 main-belt asteroids (MBAs), 4 
Jupiter Trojans, 5 Centaurs, and 7 trans-Neptunian objects 
(TNOs). 

These classes of objects present very different physical 
properties such as sizes, regolith grain size, average value of 
the thermal inertia, and composition. Other important dif¬ 
ferences are their average surface temperature due to their 
very different heliocentric distances and orbital elements. 
The illumination and observation geometry are also diverse 
for different classes of objects. Eor instance, for TNOs and 
MBAs the phase angle of observation from Earth and Earth¬ 
like orbits is typically between a few and a few tens of de¬ 
grees, respectively. On the other hand, NEAs can be ob¬ 
served under a much wider range of phase angles than can 
approach hundred degrees and more Muller (see also 2002). 
A special care should be used in these cases to explicitly 
calculate the heat diffusion in craters instead of using the 
approximation of Eq. 9. 

6.1. Ground truth from space missions to asteroids 

A number of asteroids have been, or will be, visited 
by spacecraft, providing ground-truth for the application of 
TPMs to remote-sensing thermal-infrared data. 

(21) Lutetia: based on ground-based data and a TPM, 
Mueller et ak (2006) measured Lutetia’s effective diame¬ 
ter and pv to within a few percent of the later Rosetta re¬ 
sult (Sierks et ak, 2011). Their thermal-inertia constraint 
(r < 50 J m'^s’^^^K'^) was refined by O’Rourke et ak 
(2012) based on the Rosetta shape model and more than 
70 thermal-infrared observations obtained from the ground, 
Spitzer, Akari, and Herschel: E = 5 J m‘^s‘'^^K‘'with a 
high degree of surface roughness. Keihm et ak (2012) used 
MIRO aboard Rosetta to obtain a surface thermal inertia 
< 30 J m'^s’^^^K'^. The low thermal inertia can be ex¬ 
plained by a surface covered in fine regolith; Gundlach and 
Blum (2013) infer a regolith grain size of about 200 pm. 
The study of the morphology of craters by Vincent et ak 
(2012) indicates abundant, thick (600 m), and very fine re- 
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Table 2 

PUBLISED THERMAL INERTIA VALUES 


Number 

Name 

D 

(km) 

Ao 

(km) 

r 

(SI) 

Ar 

(SI) 

Tax 

r 

(au) 

Ref. 

Number 

Name 

D 

(km) 

Ao 

(km) 

r 

(SI) 

Ar 

(SI) 

Tax 

r 

(au) 

Ref. 

1 

Ceres 

923 

20 

10 

10 

C 

2.767 

1 

1620 

Geographos 

5.04 

0.07 

340 

120 

S 

1.1 

12 

2 

Pallas 

544 

43 

10 

10 

B 

2.772 

1 

1862 

Apollo 

1.55 

0.07 

140 

100 

Q 

1.0 

13 

3 

Juno 

234 

11 

5 

5 

S 

2.671 

1 

2060 

Chiron 

142 

10 

4 

4 

B/Cb 

8-15 

14 

4 

Vesta 

525 

1 

20 

15 

V 

2.3 

2 

2060 

Chiron 

218 

20 

5 

5 

B/Cb 

13 

15 

16 

Psyche 

244 

25 

125 

40 

M 

2.7 

3 

2363 

Cebriones 

82 

5 

7 

7 

D 

5.2 

16 

21 

Lutetia 

96 

1 

5 

5 

M 

2.8 

4 

2867 

Steins 

4.92 

0.4 

150 

60 

E 

2.1 

17 

22 

Kalliope 

167 

17 

125 

125 

M 

2.3 

5 

2867 

Steins 

5.2 

1 

210 

30 

E 

2.1 

18 

32 

Pomona 

85 

1 

70 

50 

S 

2.8 

6 

8405 

Asbolus 

66 

4 

5 

5 

- 

7.9 

19 

41 

Daphne 

202 

7 

25 

25 

Ch 

2.1 

7 

3063 

Makhaon 

116 

4 

15 

15 

D 

4.7 

16 

44 

Nysa 

81 

1 

120 

40 

E 

2.5 

6 

10199 

Chariklo 

236 

12 

1 

1 

D 

13 

14 

45 

Eugenia 

198 

20 

45 

45 

C 

2.6 

5 

10199 

Chariklo 

248 

18 

16 

14 

D 

13 

13 

87 

Sylvia 

300 

30 

70 

60 

P 

2.7 

5 

25143 

Itokawa 

0.32 

0.03 

700 

100 

S 

1.1 

8 

107 

Camilla 

245 

25 

25 

10 

P 

3.2 

5 

25143 

Itokawa 

0.320 

0.029 

700 

200 

S 

1.1 

20 

no 

Lydia 

93.5 

3.5 

135 

65 

M 

2.9 

6 

29075 

1950 DA 

1.30 

0.13 

24 

20 

M 

1.7 

21 

115 

Thyra 

92 

2 

62 

38 

S 

2.5 

6 

33342 

1998 WT 24 

0.35 

0.04 

200 

100 

E 

1.0 

8 

121 

Hermione 

220 

22 

30 

25 

Ch 

2.9 

5 

50000 

Quaoar 

1082 

67 

6 

4 

- 

43 

15 

130 

Elektra 

197 

20 

30 

30 

Ch 

2.9 

5 

54509 

YORP 

0.092 

0.010 

700 

500 

S 

1.1 

8 

277 

Elvira 

38 

2 

250 

150 

S 

2.6 

6 

55565 

2002 AWigy 

700 

50 

10 

10 

- 

47 

22 

283 

Emma 

135 

14 

105 

100 

P 

2.6 

5 

90377 

Sedna 

995 

80 

0.1 

0.1 

- 

87 

23 

306 

Unitas 

56 

1 

180 

80 

S 

2.2 

6 

90482 

Orcus 

968 

63 

1 

1 

- 

48 

13 

382 

Dodona 

75 

1 

80 

65 

M 

2.6 

6 

99942 

Apophis 

0.375 

0.014 

600 

300 

Sq 

1.05 

24 

433 

Eros 

17.8 

1 

150 

50 

S 

1.6 

8 

101955 

Bennu 

0.495 

0.015 

650 

300 

B 

1.1 

25 

532 

Herculina 

203 

14 

10 

10 

S 

2.772 

1 

101955 

Bennu 

0.49 

0.02 

310 

70 

B 

1.1 

26 

617 

Patroclus 

106 

11 

20 

15 

P 

5.9 

9 

136108 

Haumea 

1240 

70 

0.3 

0.2 

- 

51 

27 

694 

Ekard 

109.5 

1.5 

120 

20 

- 

1.8 

6 

162173 

1999 JU 3 

0.87 

0.03 

400 

200 

C 

1.4 

28 

720 

Bohlinia 

41 

1 

135 

65 

s 

2.9 

6 

175706 

1996 FG 3 

1.71 

0.07 

120 

50 

c 

1.4 

29 

956 

Elisa 

10.4 

0.8 

90 

60 

- 

1.8 

10 

208996 

2003 AZ84 

480 

20 

1.2 

0.6 

- 

45 

27 

1173 

Anchises 

136 

15 

50 

20 

p 

5.0 

11 

308635 

2005 YU 55 

0.306 

0.006 

575 

225 

c 

1.0 

30 

1580 

Betulia 

4.57 

0.46 

180 

50 

c 

1.1 

8 

341843 

2008 EV 5 

0.370 

0.006 

450 

60 

c 

1.0 

31 


References: [1] Muller and Lagerros (1998), [2] Leyrat et al. (2012), [3] Matter et al. (2013), [4] O’Rourke et al. (2012), [5] Marchis et al. (2012), [6] Delbo and Tanga 
(2009), [7] Matter et al. (2011), [8] Mueller (2007), [9] Mueller et al. (2010), [10] Lim et al. (2011), [11] Homer et al. (2012), [12] Rozitis and Green (2014), [13] Rozitis 
et al. (2013), [14] Groussin et al. (2004), [15] Fornasier et al. (2013), [16] Fernandez et al. (2003), [17] Lamy et al. (2008), [18] Leyrat et al. (2011), [19] Fernandez et al. 
(2002), [20] Muller et al. (2014a), [21] Rozitis et al. (2014), [22] Cruikshank et al. (2005), [23] Pal et al. (2012), [24] Muller et al. (2014b), [25] Muller et al. (2012), [26] 
Emery et al. (2014), [27] Lellouch et al. (2013), [28] Muller et al. (2011), [29] Wolters et al. (2011), [30] Muller et al. (2013), [31] Ali-Lagoa et al. (2014), Note: for Ceres, 

Pallas, Juno, and Herculina r is assumed equal to the semimajor asxis of the orbit. 
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golith, confirming the TPM results. 

(433) Eros was studied by the NASA NEAR-Shoemaker 
space mission that allowed determination of the shape and 
size of this asteroid (mean radius of 8.46 km with a mean 
error of 16 m; Thomas et ah, 2002). Mueller (2007) per¬ 
formed a TPM analysis of the ground-based thermal in¬ 
frared data by Harris and Davies (1999), obtaining a best-fit 
diameter of 17.8 km that is within 5% of the Thomas et al. 
(2002) result of 16.9 km, and T in the range 100 - 200 
J m‘^s‘^^^K '. The latter value, in agreement with TPM re¬ 
sults of Lebofsky and Rieke (1979), implies coarser surface 
regolith than that on the Moon and larger asteroids (see, 
e.g., Mueller, 2007; Delbo et al., 2007). From the value 
of r of Mueller (2007), Gundlach and Blum (2013) calcu¬ 
lated a 1-3 mm typical regolith grain size for Eros. Optical 
images of the NEAR-Schoemaker landing site at a resolu¬ 
tion of about 1 cm/pixel (Veverka et al., 2001) show very 
smooth areas at the scale of the camera spatial resolution 
(Fig. 8), likely implying mm or sub-mm grain size regolith, 
consistent with TPM results. 

(25143) Itokawa physical properties were derived in- 
situ by the JAXA sample-return mission Hayabusa, allow¬ 
ing us to compare the size, albedo and regolith nature de¬ 
rived from the TPMs with spacecraft results. Muller et al. 
(2014a) show an agreement within 2% between the size and 
the geometric visible albedo inferred from TPM analysis 
of thermal-infrared data and the value of the corresponding 
parameters from Hayabusa data. The TPM thermal inertia 
value for Itokawa is around 750 J m‘^s‘*^^K *, significantly 
higher than the value of our Moon (about 50 J m'^s’^^^K'^) 
and of other large main belt asteroids including (21) Lutetia, 
implying a coarser regolith on this small NEA. The corre¬ 
sponding average regolith grain size according to Gundlach 
and Blum (2013) is ~ 2 cm. Hayabusa observations from 
the optical navigation camera (ONC-T), obtained during the 
descent of the spacecraft to the “Muses Sea” region of the 
asteroid, reveal similar grain sizes, at a spatial resolution of 
up to 6 mm/pixel. In particular, Yano et al. (2006) describe 
“Muses Sea” as composed of numerous size-sorted granular 
materials ranging from several centimeters to subcentime¬ 
ter scales. Itokawa’s regolith material can be classified as 
’’gravel”, larger than submillimeter regolith powders filling 
in ponds on (433) Eros (Fig. 8). 

It is worth pointing out, however, that “Muses Sea” is 
not representative of Itokawa’s surface as a whole. Rather, 
it was selected as a touchdown site because, in earlier 
Hayabusa imaging, it appeared as particularly smooth (min¬ 
imizing operational danger for the spacecraft upon touch¬ 
down) and apparently regolith rich (maximizing the chance 
of sampling regolith). Grain sizes measured at “Muses Sea” 
are therefore lower limits on typical grain sizes rather than 
values typical for the surface as a whole. 

6.2. Thermal inertia of large and small asteroids 

An inverse correlation between E and D was noticed by 
Delbo et al. (2007), then updated Delbo and Tanga (2009) 



Fig. 8.— Higher P-values correspond to coarser regoliths. (A) 
Close-up image of (433) Eros from the NASA NEAR Shoemaker 
mission reveals coarse regolith with grain size in the mm-range 
(adapted from Veverka et al., 2001). The value of T is ~150 
J m"^s‘*^^K'* for Eros. (B) Image from the JAXA Hayabusa mis¬ 
sion (from Yano et al., 2006) of the surface of (25143) Itokawa dis¬ 
playing gravel-like regolith. The value of F is ~750 J m"^s‘*^^K * 
for Itokawa. 

and Capria et al. (2014). This supported the intuitive view 
that large asteroids have, over many hundreds of millions of 
years, developed substantial insulating regolith layers, re¬ 
sponsible for the low values of their surface thermal inertia. 
On the other hand, much smaller bodies, with shorter col- 
lisional lifetimes (Marchi et al., 2006; Bottke et al., 2005, 
and references therein), have less regolith, and or larger re¬ 
golith grains (less mature regolith), and therefore display a 
larger thermal inertia. 

In the light of the recently published values of E (Tab. 2), 
said inversion correlation between E and D is less clear, in 
particular, when the values of the thermal inertia are tem¬ 
perature corrected (Fig. 9). However, the E vs U distribu¬ 
tion of ZJ > 100 km (large) asteroids is different than that 
of ZJ < 100 km (small) asteroids. Small asteroids typically 
have higher E-values than large asteroids, which present a 
large scatter of E-values, ranging from a few to a few hun¬ 
dreds J m'^s '^^K '. This is a clear indication of a diverse 
regolith nature amongst these large bodies. A shortage of 
low E values for small asteroids is also clear, with the no¬ 
table exception of 1950 DA, which has an anomalously low 
E-value compared to other NEAs of similar size (Rozitis 
et al., 2014). 

Fig. 9 also shows previously unnoticed high-thermal- 
inertia C types, maybe related to CR carbonaceous chon¬ 
drites, which contain abundant metal phases. We also 
note that all E types in our sample appear to have a size- 
independent thermal inertia. 

6.3. Very low E-values 

We also note that the some of the C-complex outer 
main-belt asteroids and Jupiter Trojans have very low ther¬ 
mal inertia in the range between a few and a few tens of 
J m'^s '^^K"'. In order to reduce the thermal inertia of a ma¬ 
terial by at least one order of magnitude (from the lowest 
measured thermal inertia of a meteorite, ~650 J m‘^s‘'^^K ' 
at 200 K (Opeil et al., 2010), to the typical values for these 
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Fig. 9.— r values vs. D from Tab. 2 for different taxonomic 
types (see key). Top plot: original measurements, bottom plot: T 
corrected to 1 au heliocentric distance for temperature dependent 
thermal inertia assuming Eq. 13 and the heliocentric distance at the 
time of thermal infrared observations reported in Tab. 2. Trojans, 
Centaurs and trans-Neptunian objects are not displayed. 

large asteroids (Tab. 2 and Fig. 9), a very large porosity 
(>90%) of the first few mm of the regolith is required 
(Vemazza et al., 2012). This is consistent with the dis¬ 
covery that emission features in the mid-infrared domain 
(7-25 pm, Fig. 4) are rather universal among large aster¬ 
oids and Jupiter Trojans (Vernazza et al., 2012), and that 
said features can be reproduced in the laboratory by sus¬ 
pending meteorite and/or mineral powder (with grain sizes 
< 30 pm) in IR-transparent KBr (potassium bromide) pow¬ 
der (Vernazza et al., 2012). As KBr is not supposed to 
be present on the surfaces of these minor bodies, regolith 
grains must be ’’suspended” in void space likely due to co¬ 
hesive forces and/or dust levitation. On the other hand, 
radar data indicate a significant porosity (40-50 %) of the 
first ^1 m of regolith (Magri et al., 2001; Vernazza et al., 
2012 ), indicating decreasing porosity with increasing depth 
(see Fig. 5 of Vemazza et al., 2012, for a regolith schemat¬ 
ics). 


6.4. Average thermal inertia of asteroid populations 

As described before, the thermal inertia of an aster¬ 
oid can be directly derived by comparing measurements 
of its thermal-infrared emission to model fluxes generated 
by means of a TPM. Typically, more than one observation 
epoch is required to derive the thermal inertia, in order to 
’’see” the thermal emission from different parts of the as¬ 
teroid’s diurnal temperature distribution. Unfortunately, the 
large majority of minor bodies for which we have thermal- 
infrared observations have been observed at a single epoch 
and/or information about their gross shape and pole orienta¬ 
tion is not available, precluding the use of TPMs. However, 
if one assumes the thermal inertia to be roughly constant 
within a population of asteroids (e.g., NBAs) one can use 
observations of different asteroids under non-identical illu¬ 
mination and viewing geometries, as if they were from a 
unique object. Delbo et al. (2003) noted that qualitative in¬ 
formation about the average thermal properties of a sample 
of NBAs could be obtained from the distribution of the re¬ 
values of the sample as a function of the phase angle, a. 
Delbo et al. (2007) and Bellouch et al. (2013) developed 
a rigorous statistical inversion method, based on the com¬ 
parison of the distributions of published NBATM p-values 
vs a, or vs. r with that of a synthetic population of aster¬ 
oids generated through a TPM, using realistic distributions 
of the input TPM parameters such as the rotation period, 
the aspect angle etc. Delbo et al. (2007) found that the av¬ 
erage thermal inertia value for km-sized NBAs is around 
200 J m'^s '^^K'^. The average thermal inertia of binary 
NBAs is higher than that of non-binary NBAs, possibly in¬ 
dicating a regolith-depriving mechanism for the formation 
of these bodies (Delbo et al., 2011). The same authors also 
found that NBAs with slow rotational periods {P >10 h) 
have higher-than-average thermal inertia. Brom a sample 
of 85 Centaurs and trans-Neptunian objects observed with 
Spitzer/MIPS and HerscheFPACS, Bellouch et al. (2013) 
found that surface roughness is significant, a mean ther¬ 
mal inertia F = 2.5 ± 0.5 J m'^s '^^K ^ and a trend to¬ 
ward decreasing F with increasing heliocentric distance. 
The thermal inertias derived by Bellouch et al. (2013) are 
2-3 orders of magnitude lower than expected for compact 
ices, and generally lower than on Saturn’s satellites or in 
the Pluto/Charon system. These results are suggestive of 
highly porous surfaces. 

6.5. Relevant astronomical and laboratory data 

Physical interpretations of thermal-inertia estimates de¬ 
pend strongly on laboratory and ground-truth measure¬ 
ments of relevant material properties. While in the Asteroid 
III era, we based interpretation of thermal inertia on Barth 
analog materials, in the last few years laboratory measure¬ 
ments were performed on asteroid analog materials, i.e., 
meteorites. Meteorite grain densities range from ~2800 
kg m'^ for CM carbonaceous chondrites to ~3700 kg m‘^ 
for enstatite chondrites (Consolmagno et al., 2006; Macke 
et al., 2010, 201 la,b). Heat capacities have been measured 
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for a wide sampling of meteorites by Consolmagno et al. 
(2013), who find that values for stony meteorites are be¬ 
tween 450 and 550 J kg'^ K *, whereas C for irons tends 
to be smaller (330 - 380 J kg ^ K '). Opeil et al. (2012, 
2010 ) present thermal conductivity measurements of stony 
meteorites, finding values of 0.5 W K"' m'^ for the car¬ 
bonaceous chondrite Cold Bokkeveld to 5.5 W K"' m * for 
the enstatite chondrite Pillistfer. Their one iron meteorite 
sample has a k of 22.4 W K"' m '. They also find a linear 
correlation between and the inverse of the porosity, from 
which Opeil et al. (2012) conclude that the measured k of 
the samples is controlled more by micro-fractures than by 
composition. 

Grain size and packing, more than compositional het¬ 
erogeneity, are responsible for different thermal inertias of 
different surfaces. This also explains why TPMs are capa¬ 
ble of deriving asteroid physical parameters independently 
of the asteroid mineralogy. Conduction between grains is 
limited by the area of the grain contact (Piqueux and Chris¬ 
tensen, 2009b,a). As grain size decreases to diameters less 
than about a thermal skin depth (few cm on most aster¬ 
oids), conduction is more and more limited (e.g., Presley 
and Christensen, 1997b). On bodies with atmospheres, con¬ 
duction through the air in pores can often efficiently trans¬ 
port heat. On airless bodies, however, radiation between 
grains, which is not very efficient, particularly at low T 
(e.g., Gundlach and Blum, 2012), is the only alternative to 
conduction across contacts (Fig. 10). Considering these two 
modes of energy transport and their dependence on grain 
size, Gundlach and Blum (2013) developed an analytical 
approach for determining grain size from thermal inertia 
measurements. They incorporated the measurements of ma¬ 
terial properties of meteorites measured above along with 
results of their own laboratory of heat transport in dusty 
layers. Additional laboratory measurements of conductiv¬ 
ities of powdered meteorites under high vacuum would be 
valuable for more precise interpretation of asteroid thermal 
inertias. 

The classic opportunity for ground-truth thermal mea¬ 
surements came with the Apollo missions. Astronauts 
on Apollo 15 and 17 carried out bore-hole style temper¬ 
ature measurements to depths of 1.4 m below the surface 
on Apollo 15 and 2.3 m below the surface on Apollo 17 
(Langeseth and Keihm, 1977; Vaniman et al., 1991). Ther¬ 
mal conductivity of about 0.001 W K"' m ' was found in the 
top 2 to 3 cm of the lunar regolith, increasing to about 0.01 
W K'^ m * over the next few cm, then to values as high as 
2 W K * m * deeper into the surface where the regolith ap¬ 
pears to have been very compacted (Langeseth and Keihm, 
1977). Low thermal inertias derived from remote thermal 
infrared measurements (e.g. Wesselink, 1948b; Vasavada 
et al., 2012) agree with the very low k in the topmost few cm 
of the lunar surface, and the Apollo measurements provide 
the necessary ground-trath for interpreting such low ther¬ 
mal inertias as very fine-grained, ’’fluffy” regolith. These 
measurements fostered, for instance, development of de¬ 
tailed models of lunar regolith (Keihm, 1984). Detailed 



Bottom of the Dust Layer 


—► : Heat Transport within the Solid Material 
—*■ : Heat Transport due to Radiation 
—*■ : Heat Transport due to Gas Diffusion 
^Suif ' ^t‘tf‘nte Temperature of the Dust Layer 
^Bottom' Temperature of the Dust Layer 

Fig. 10.— Diagram of the modes of heat transport in regoliths. 
On airless bodies, heat can flow by conduction through grain 
boundaries (solid line) or by radiation between grains (dashed 
line). The dotted line showing transport by gas diffusion is not 
relevant to asteroid surfaces. From Gundlach and Blum (2012). 

thermal infrared observations and thermal models of the 
lunar regolith allows today estimating the subsurface rock 
abundance (e.g., Bandfield et al., 2011), allowing geologi¬ 
cal studies of the regolith production rate. 

6.6. Dependence of F with depth 

The depth dependence of typical asteroid regolith prop¬ 
erties is poorly constrained at this point, which is why 
physical constants are typically assumed to be constant 
with depth. MIRO observations of (21) Lutetia, how¬ 
ever, showed the existence of a top layer with F < 30 
J m'^s '^^K *, while the thermal inertia of subsurface mate¬ 
rial appears to increase with depth much like on the Moon 
(Keihm et al., 2012). 

6.7. Infrared linih brightening 

Recent modeling and observations show that, contrary 
to expectation, the flux enhancement measured in disk- 
integrated observations of the sunlit side of an asteroid (e.g., 
Lebofsky et al., 1986) is dominated by limb surfaces rather 
than the subsolar region (Rozitis and Green, 2011; Keihm 
et al., 2012). This suggests that for the sunlit side of an 
asteroid, sunlit surfaces directly facing the observer in situ¬ 
ations where they would not be if the surface was a smooth 
flat one are more important than mutual self-heating be¬ 
tween interfacing facets raising their temperatures. Figure 9 
of Rozitis and Green (2011) pictures this effect for a Gaus¬ 
sian random surface during sunrise viewed from different 
directions. The thermal flux observed is enhanced when 
viewing hot sunlit surfaces (i.e.. Sun behind the observer), 
and is reduced when viewing cold shadowed surfaces (i.e.. 
Sun in front of the observer). 

Jakosky et al. (1990) also studied the directional thermal 
emission of Earth-based lava flows exhibiting macroscopic 
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roughness. They found that enhancements in thermal emis¬ 
sion were caused by viewing hot sunlit sides of rocks and 
reductions were caused by viewing cold shadowed sides of 
rocks. This agrees precisely with the model and adds fur¬ 
ther evidence that thermal infrared beaming is caused by 
macroscopic roughness rather than microscopic roughness. 

The effect of limb brightening has also been measured 
from disk-resolved thermal infrared data (<5 pm) acquired 
during sunrise on the nucleus of the comet 9P/Tempel 1 by 
the Deep Impact NASA space mission (Davidsson et ak, 
2013), and from VIRTIS and MIRO measurements of the 
asteroid (21) Lutetia (Keihm et ak, 2012). 

6.8. Asteroid thermal inertia maps 

Disk-resolved thermal infrared observations, in the range 
between 4.5 -5.1 pm, were provided by the instrument VIR 
(De Sanctis et ak, 2012) on board of the NASA DAWN 
(Russell et ak, 2012) spacecraft (Capria et ak, 2014, and 
references therein). Form TPM analysis of VIR measure¬ 
ments, Capria et ak (2014) obtained a map of the roughness 
and the thermal inertia of Vesta. The average thermal in¬ 
ertia of Vesta is 30 ± 10 J m’^s '^^K"', which is in good 
agreement with the values found by ground-based observa¬ 
tions (Muller and Lagerros, 1998; Chamberlain et ak, 2007; 
Leyrat et ak, 2012). The best analog is probably the sur¬ 
face of the Moon, as depicted by Vasavada et ak (2012) and 
Bandfield et ak (2011): a surface whose thermal response is 
determined by a widespread layer of dust and regolith with 
different grain sizes and density increasing toward the in¬ 
terior. Exposed rocks are probably scarce or even absent. 
Capria et ak (2014) also show that Vesta cannot be consid¬ 
ered uniform from the point of view of thermal properties. 
In particular, they found that the thermal inertia spatial dis¬ 
tribution follows the global surface exposure age distribu¬ 
tion, as determined by crater counting in Raymond et ak 
(2011), with higher thermal inertia displayed by younger 
terrains and lower thermal inertia in older soils. 

Capria et ak (2014) also found higher-than-average ther¬ 
mal inertia terrain units located in low-albedo regions that 
contain highest abundance of OH, as determined by the 2.8 
pm band depth (De Sanctis et ak, 2012). These terrains are 
associated with the dark material, thought to be delivered by 
carbonaceous chondrite like asteroids that have impacted 
Vesta at low velocity. Note that in general (carbonaceous 
chondrites) have lower densities and lower thermal conduc¬ 
tivity (Opeil et ak, 2010) than basaltic material, which con¬ 
stitute the average Vestan terrain. This consideration would 
point to a lower thermal inertia rather than a higher one, as 
observed on Vesta. Capria et ak (2014) conclude that the 
factor controlling the thermal inertia in these areas could be 
the degree of compaction of the uppermost surface layers, 
which is higher than in other parts of the surface. 

6.9. Thermal inertia of metal-rich regoliths 

In principle, the composition of the regolith and not 
only its average grain size and the degree of compaction 


also affects the thermal inertia of the soil (Gundlach and 
Blum, 2013). For instance iron meteorites have a higher 
thermal conductivity than ordinary and carbonaceous chon¬ 
drites (Opeil et ak, 2010). We thus expect that a metal iron 
rich regolith displays a higher thermal inertia than a soil 
poor of this component. Harris and Drube (2014) compared 
values of the NEATM p-parameter derived from WISE data 
with asteroid taxonomic classifications and radar data, and 
showed that the p-value appears to be a useful indicator of 
asteroids containing metal. Matter et ak (2013) performed 
interferometric observations with MIDI of the ESO-VLTI 
in thermal infrared of (16) Psyche and showed that Psy¬ 
che has a low surface roughness and a thermal inertia value 
around 120 ± 40 J m'^s '^^K ', which is one of the higher 
values for an asteroid of the size of Psyche (~ 200 km). 
This higher than average thermal inertia supports the evi¬ 
dence of a metal-rich surface for this body. 

7. EFFECTS OF TEMPERATURES ON THE SUR¬ 
FACE OF ASTEROIDS 

7.1. Thermal cracking 

The surface temperature of asteroids follows a diur¬ 
nal cycle (see Fig. 2) with typically dramatic tempera¬ 
ture changes as the Sun rises or sets. The resulting, re¬ 
peated thermal stress can produce cumulative damage on 
surface material due to opening and extension of micro¬ 
scopic cracks. This phenomenon is known as thermal fa¬ 
tigue (Delbo et ak, 2014). 

Growing cracks can lead to rock break-up when the num¬ 
ber of temperature cycles is large enough. For typical aster¬ 
oid properties, this process is a very effective mechanism 
for comminuting rocks and to form fresh regolith (Delbo 
et ak, 2014). For cm-sized rocks on an asteroid 1 au from 
the Sun, thermal fragmentation is at least an order of magni¬ 
tude faster than comminution by micrometeoroid impacts, 
the only regolith-production mechanism previously consid¬ 
ered relevant (Horz and Cintala, 1997; Hoerz et ak, 1975). 

The efficiency of thermal fragmentation is dominated by 
the amplitude of the temperature cycles and by the temper¬ 
ature change rate (Hall and Andre, 2001), which in turn de¬ 
pend on heliocentric distance, rotation period, and the sur¬ 
face thermal inertia. The rate of thermal fragmentation in¬ 
creases with decreasing perihelion distance: at 0.14 au from 
the Sun, thermal fragmentation may erode asteroids such as 
(3200) Phaethon and produce the Geminids (Jewitt and Li, 
2010), whereas in the outer Main Belt this process might 
be irrelevant. Thermal fragmentation of surface boulders is 
claimed by Dombard et ak (2010) to be source of fine re¬ 
golith in the so-called ’’ponds” on the asteroid (433) Eros. 
Production of fresh regolith originating in thermal fatigue 
fragmentation may be an important process for the rejuve¬ 
nation of the surfaces of near-Earth asteroids (Delbo et ak, 
2014). 

Thermal cracking is reported on other bodies, too: on 
Earth, particularly in super-arid environments (Hall, 1999; 
Hall and Andre, 2001), on the Moon (Levi, 1973; Duen- 
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nebier and Sutton, 1974), Mercury (Molaro and Byrne, 
2012), Mars (Viles et al., 2010), and on meteorites (Levi, 
1973). Moreover, Tambovtseva and Shestakova (1999) sug¬ 
gest that thermal cracking could be an important process in 
the fragmentation and splitting of kilometer-sized comets 
while in the inner solar system. Furthermore, Capek and 
Vokrouhlicky (2010) initially proposed that slowly rotating 
meteoroids or meteoroids that have spin vector pointing to¬ 
wards the sun can be broken up by thermal cracking. In 
a further development of their model, Capek and Vokrouh¬ 
licky (2012) showed that as the meteoroid approaches the 
Sun, the stresses first exceed the material strength at the sur¬ 
face and create a fractured layer. If inter-molecular forces 
(e.g., Rozitis et al., 2014) are able to retain the surface layer, 
despite the competing effects of thermal lifting and centrifu¬ 
gal forces, the particulate surface layer is able to thermally 
shield the core, preventing any further damage by thermal 
stresses. 

7.2. Sun-driven heating of near-Earth asteroids and 
meteoroids 

It is known that heating processes can affect the physical 
properties of asteroids and their fragments, the meteorites 
(see, e.g., Keil, 2000). 

Internal heating due to the decay of short-lived radionu¬ 
clides was considered early on (Grimm and McSween, 
1993). March! et al. (2009) discuss close approaches to the 
Sun as an additional surface-altering heating mechanism. 
In the present near-Earth asteroid population, the fraction 
of bodies with relatively small perihelion (q) is very small: 
about 1/2, 1/10, and 1/100 of the population of currently 
known near-Earth objects (11,000 as of the time of writ¬ 
ing) have a perihelion distance below 1, 0.5, and 0.25 au, 
where maximum temperature are exceeding 400, 550, and 
780 K, respectively (see Eig. 11). However, dynamical sim¬ 
ulations show that a much larger fraction of asteroids had 
small perihelion distances for some time, hence experienc¬ 
ing episodes of strong heating in their past (March! et al., 
2009). Eor instance, the asteroid 2004 LG was approaching 
the Sun to within only ~5.6 solar radii some 3 ky ago, and 
its surface was baked at temperatures of 2500 K (Vokrouh¬ 
licky and Nesvorny, 2012). 

Solar heating has a penetration depth of typically a few 
cm (see Eq. 5 and Spencer et al., 1989). Organic compo¬ 
nents found on meteorites break up at temperatures as low 
as 300-670 K (see Eig. 11 and Kebukawa et al., 2010; Erost 
et al., 2000; Huang et al., 1994), thus solar heating can re¬ 
move these components from asteroid surfaces. 

7.3. Thermal metamorphism of meteorites 

Radiative heating from the Sun has been invoked as a 
mechanism for the thermal metamorphism of metamorphic 
CK carbonaceous chondrites (Chaumard et al., 2012). The 
matrix of these chondrites shows textures consistent with a 
transient thermal event during which temperatures rose be¬ 
tween 550 and 950 K. The inferred duration of these events 



Perihelion Distance (q; AU) 

Eig. 11.— Surface temperature of an asteroid or meteoroid as a 
function of the distance from the Sun. Vertical arrows indicate the 
threshold temperature for the thermal alteration/desiccation for a 
variety of chemical compounds discussed in the text (see Delbo 
and Michel, 2011, and references therein for further information). 
The temperature range for thermal metamorphism of the CK chon¬ 
drites is from Chaumard et al. (2012). 

is of the order of days to years, much longer than the time 
scale of shock events but shorter than the time scale for 
heating by the decay of radiogenic species such as ^®A1 
(e.g., Kallemeyn et al., 1991). 

7.4. Subsurface ice sublimation 

Observational evidence for the presence of ice on aster¬ 
oid surfaces stems from the discovery of main belt comets 
(MBCs; Hsieh and Jewitt, 2006), the localized release of 
water vapor from the surface of (1) Ceres (Kiippers et al., 
2014), and the detection of spectroscopic signatures inter¬ 
preted as water ice frost on the surface of (24) Themis 
(Rivkin and Emery, 2010; Campins et al., 2010) and of (65) 
Cybele (Licandro et al., 2011). 

The lifetime of ices on the surface and in the subsurface 
depends strongly on temperature. TPMs have been used to 
estimate these temperatures. This requires a modification 
of the “classical” TPM as presented in section 4, such that 
heat conduction is coupled with gas diffusion (Schorghofer, 
2008; Capria et al., 2012; Prialnik and Rosenberg, 2009). 
The referenced models assume a spherical shape. As for the 
interior structure, Capria et al. (2012); Prialnik and Rosen¬ 
berg (2009) assume a comet-like structure, i.e., an inti¬ 
mate mixture of ice and dust throughout the entire body, 
while Schorghofer (2008) consider an ice layer underneath 
a rocky regolith cover. Sublimation of ice and the trans¬ 
port of water molecules through the fine-grained regolith is 
modeled in all cases. 

All authors agree that water ice exposed on asteroid sur¬ 
faces sublimates completely on timescales much shorter 
than the age of the Solar System. Therefore, asteroid sur¬ 
faces were expected to be devoid of water ice, contrary to 
the observational evidence quoted above. However, water 
ice can be stable over 4.5 Gy in the shallow subsurface, at a 
depth of ~1-10 m. In particular, Eanale and Salvail (1989) 
showed that ice could have survived in the subsurface at the 
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polar regions of Ceres. Large heliocentric distances, slow 
rotation, and a fine-grained regolith leading to low thermal 
conductivity and short molecular free path, all favor the sta¬ 
bilization of subsurface water ice (Schorghofer, 2008). The 
same authors conclude that rocky surfaces, in contrast to 
dusty surfaces, are rarely able to retain ice in the shallow 
subsurface. 

To be observable on the surface, buried ice most be ex¬ 
posed. Campins et al. (2010) describe several plausible 
mechanisms such as impacts, recent change in the obliquity 
of the spin pole, and daily or orbital thermal pulses reaching 
a subsurface ice layer. 

8. FUTURE CHALLENGES FOR TPMs 

The Spitzer and WISE telescopes have opened a new era 
of asteroid thermal-infrared observations and the exploita¬ 
tion of their data through TPMs has just begun (e.g., Ali- 
Lagoa et al., 2014; Rozitis et al., 2014; Emery et al., 2014). 
At the moment, the limiting factor is the availability of ac¬ 
curate asteroid shape models. However, optical-wave length 
all-sky surveys such as PanSTARRS, LSST, and Gaia are 
expected to produce enormous photometric data sets lead¬ 
ing to thousands of asteroid models. We envision the avail¬ 
ability of thousands of thermal-inertia values in some years 
from now, enabling more statistically robust studies of ther¬ 
mal inertia as a function of asteroid size, spectral class, 
albedo, rotation period, etc. 

For instance, the distribution of T within asteroid fam¬ 
ilies will be crucial in the search of evidence of aster¬ 
oid differentiation; asteroid formation models and mete¬ 
orite studies suggest that hundreds of planetesimals experi¬ 
enced complete or partial differentiation. An asteroid fam¬ 
ily formed from the catastrophic disruption of such a dif¬ 
ferentiated asteroid should contain members corresponding 
to the crust, the mantel and the iron core. However, the 
observed spectra and albedos are very homogeneous across 
asteroid families. Thermal inertia might help in separating 
iron rich from iron poor family members, supposedly orig¬ 
inating respectively from the core and mantle of the differ¬ 
entiated parent body (e.g.. Matter et al., 2013; Harris and 
Drube, 2014). 

At any size range. Fig. 9 shows an almost tenfold vari¬ 
ability in thermal inertia, corresponding to difference in av¬ 
erage regolith grain size of almost two orders of magnitude 
Gundlach and Blum (2013). For small near-Earth aster¬ 
oids, this could be due to a combination of thermal cracking 
(Delbo et al., 2014), regolith motion (Murdoch et al., 2015), 
and cohesive forces (Rozitis et al., 2014). Faster rotation 
periods allow more thermal cycles, which then enhances 
thermal fracturing. It also encourages regolith to move to¬ 
wards the equator where the gravitational potential is at its 
lowest (Walsh et al., 2008). And for the extremely fast rota¬ 
tors, large boulders/rocks could be selectively lofted away, 
because they stick less well to the surface than smaller par¬ 
ticles. For D > 100 km sized asteroids, F-values might be 
help to distinguish between primordial and more recently 


re-accumulated asteroids. The former had ~ 4 Gy of re¬ 
golith evolution, the latter have a less developed and there¬ 
fore coarser regolith. 

The high-precision thermal-infrared data of WISE and 
Spitzer pose new challenges to TPMs, as model uncertain¬ 
ties are now comparable to the uncertainty of the measured 
flux. This will become even more important with the launch 
of the James Webb Space Telescope (JWST). In particular, 
the accuracy of the shape models might represent a limiting 
factor (e.g., Rozitis and Green, 2014). The next challenge 
will be to allow the TPM to optimize the asteroid shape. 
This seems to be possible as the infrared photometry is also 
sensitive to shape, provided good-quality thermal data are 
available (Durech et al., 2015). 

New interferometric facilities, such as MATISSE, LBTI, 
and ALMA, will become available in the next years requir¬ 
ing TPMs to calculate precise disk-resolved thermal fluxes 
(Durech et al., 2015). The wavelengths of ALMA, simi¬ 
lar to those of MIRO, will allow to measure the thermal- 
infrared radiation from the subsoil of asteroids, thus pro¬ 
viding further information about how thermal inertia varies 
with depth. 

Certainly, constraining roughness is one of the future 
challenges for TPMs. To do so from disk-integrated data re¬ 
quires a range of wavelengths and solar phase angles. Low 
phase angle measurements are enhanced by beaming whilst 
high phase angle measurements are reduced by beaming. 
In particular, shorter wavelengths are affected more than 
longer wavelengths. 

Moreover, the future availability of precise sizes and 
cross sections of asteroids from stellar occultation timing 
(Tanga and Delbo, 2007), combined with shape informa¬ 
tion derived from lightcurve inversion (Durech et al., 2015) 
will allow to remove the need to constrain the object size 
from TPM analysis. linfrared fluxes will thus be converted 
into highly reliable thermal inertia and roughness values. 
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